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In a gauge invariant formulation of the molecular electric 
dipole-photon interaction, the rigorous coupling is strictly lin- 
ear in the photon creation and photon annihilation operators. 
The linear coupling allows for a super-radiant phase tran- 
sition as in the Hepp-Lieb formulation. A previous notion 
of a quadratic-coupling "no-go theorem" for super-radiance 
is incorrect. Also incorrect is a previous assertion that the 
dipole-photon coupling has absolutely no effect on the ther- 
mal equations of state. These dubious assertions were based 
on incorrect canonical transformations which eliminated the 
electric field (and thereby eliminated the dipole-photon inter- 
action) which is neither mathematically nor physically con- 
sistent. The correct form of the canonical transformations 
are given in this work which allows for the physical reality of 
super-radiant condensed matter phases. 

PACS: 78.60.Kn, 78.60.Fi, 78.70.-g 



I. INTRODUCTION 

In an early paper on the interaction between radiation 
and molecular dipole moments [Q, Dicke formulated a 
model which was later shown to exhibit a super-radiant 
phase transition |^,^. The notion that such phase tran- 
sitions should exist in condensed matter systems has 
been investigated in a series of papers by Preparata 
and coworkers and others [0-|9|. Different workers 
have come to somewhat different conclusions concern- 
ing super-radiant phase transitions pO|-p^. Some doubt 
has been expressed ||20|-p^ concerning the physical lab- 
oratory reality of super-radiant phase transition. The 
mathematical issues are as follows: (i) It appears, at first 
glance, that quadratic terms (in photon creation and an- 
nihilation operators) enter into the model via quadratic 
terms in the vector potential A. (ii) The quadratic terms 
in the "corrected Dicke model" appear to destroy the 
super-radiant phase transition. 

Our purpose is to show that if the dipole-field inter- 
action is treated in a gauge invariant manner 25 - 2^] 
then the interaction is strictly linear in the electric field 
E. Thus, quadratic terms are absent for purely elec- 
tric dipole-photon interactions [28] . These considerations 
render likely the physical reality of condensed matter 
super-radiant phase transitions. 

To see what is involved, let us first consider a single 
molecule with an electric dipole moment /x given by 



N 



a=l 



(1) 



where qa and Va denote respectively, the charge and posi- 
tion of the a*'' particle. Suppose that the dipole moment 
interacts with a classical electric field which is uniform 
in space, but not in time 



1 /dAU) 



(2) 



The molecule will have a time dependent Hamiltonian of 
the form 



N 



where Ka{pa) = (Pa/^'^ia) is the kinetic energy of a*'* 
particle in the molecule and V{ri, ...,r7v) is the internal 
Coulomb energy of the molecule. Employing the unitary 
operator 

U{t) = e*(M-A(t)/nc) 
in the time dependent canonical transformation 

n^oiit) = uHt)H,r,oimit) ~ ^nuHt)^ (5) 

at 



yields 



N 



= ^A'a(Pa) + "^^(ri,...,rN)-/i-E(i). (6) 

For completeness of presentation, a derivation of Eq.(5) 
is given in Appendix A. 

Note: For classical electric fields of the form in Eq.(2), 
the Hamiltonians in Eqs.(3) and (6) are rigorously equiva- 
lent. Furthermore, the Hamiltonian in Eq.(3) depends on 
linear and quadratic terms in the vector potential A(t). 
On the other hand, the Hamiltonian in Eq.(6) depends 
only on the electric field E(t) and in a strictly gauge in- 
variant and linear fashion. 

In the work which follows, both the charged particles 
and the electromagnetic field will be treated using quan- 
tum mechanics. In Sec. II we discuss a single molecule in 
the presence of the quantized electromagnetic field. It is 
proved for the gauge invariant description of the dipole 
interactions that the coupling between the electromag- 
netic field and the molecule is linear in photon creation 



1 



and annihilation operators The results for M molecules 
arc discussed in Sec. III. The resulting Hamiltonian has 
four terms: (i) a sum of single molecule Hamiltonians, 
(ii) a sum of screened Coulomb interaction potentials be- 
tween neighboring molecules, (iii) the radiation field en- 
ergy and (most importantly) (iv) a strictly linear coupling 
between the dipole moments and the electric field. Dicke- 
Preparata models are defined and explored in Sec. IV. 
In the concluding Sec.V, the importance of the super- 
radiant phase transition associated with linear couplings 
will be discussed. The errors by some previous workers 
who falsely "voided" the super-radiant phase transition 
will be discussed in detail. 



II. PHOTON-DIPOLE INTERACTIONS 



leading to Eq.(lO). 

Employing the unitary operator 

S — e^^^'-^^^'^^ 

in the time-independent canonical transformation 

n = S^HS, 

yields 

n = Hrad + Hmol{A = 0) - /X • E + Wmol 



(14) 



(15) 



(16) 



where the "self energy" contribution to the molecule is 
given by 



(17) 



Consider, as in Sec. I, a single molecule having a to- 
tal number of A'' charged particles with an electric dipole 
moment /i as in Eq.(l). The dipole is considered to inter- 
act with a quantized electromagnetic field (photons) A 
and/or E which is again uniform in space. The Hamilto- 
nian for the photons is given by 



Hrad = ^ Uwkalak, 



(7) 



where 

ak',ai =5kk'- (8) 
Furthermore, for a "quantization box" of volume O, 



and 



E^ 



(9) 



(10) 



The matter Hamiltonian for a single molecule in the elec- 
tric dipole limit is then 

N 

HmoliA) = Y,K^ (pa - ^a) + y(ri, ...r;v). (11) 

a—1 

In total, the time independent Hamiltonian for the 
molecule with dipole- quantized electromagnetic field in- 
teractions reads 



H = Hrad + Hmoli-^), 



(12) 



where Hrad is defined in Eq.(7), and Hmoi{A.) is defined 
in Eqs.(9) and (11). In the quantum electrodynamic 
model, Eq.(2) is replaced by the operator equation 



E = - 



he 



[H,A] 



(13) 



The derivation of Eq.(16) is given in Appendix B. 

Introducing the boson operators bk = iak and the Her- 
mitian conjugate b\ = — iflj. yields our final Hamiltonian 



N 

n = Y^ nujkb\bk + ^Ka(pa) + V(ri, ...rjv) 

k a=l 

+ (^)l^(efc-M)(efe-M)-M-E. 



(18) 



where 



and 



^=^I^Y.^^^^'^^k+elbl), (19) 



bk',bl 



= 6, 



kk' 



(20) 



The central result of this section is the following: 

Theorem: For the electric dipolc-clcctromagnctic field 
interaction, the resulting Hamiltonian is that of free pho- 
tons plus interaction terms linear in the photon creation 
and annihilation operators. Quadratic interaction terms 
rigorously vanish. 

Proof: See Eqs.(18) and (19), which were derived solely 
on the basis of the dipole interaction. 

Let us now consider the case of M molecules, each of 
which interact with the quantized electromagnetic field 
via the electric dipole moment. 



III. MANY MOLECULES 

We consider the case of many molecules with electric 
dipole moments interacting with electromagnetic field 
modes whose minimum wavelength Xmin is large on the 
scale of the molecular size L; i.e. Xmin » L as shown 
in Fig.l below. 
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FIG. 1. Shown are many molecules of small size L when 
compared with the minimum photon wavelength Xmin- 



The Hamiltonian for J\f molecules plus the electromag- 
netic field reads 

AT 

H = Hrad + Y.H--olAA)+ Yl ^-1- (21) 

j=l i<j<i<J^ 

The radiation Hamiltonian in the Coulomb gauge 

divA{r) = (22) 

is given by 

Hrad - ^ y" + B^) d'r. (23) 
The internal j*^ molecular Hamiltonian is given by 

N 

(p(,;,) - |a(R,)) +l^(r(i,,),...,r(^;,)), (24) 

a—l 

where Rj is the central position of the j*'* molecule. 

The internal Hamiltonian of each molecule contains the 
internal Coulomb potentials. In the last term on the 
right hand side of Eq.(21), Vij describes the Coulomb 
interactions between two different molecules. The long 
ranged part of this potential is of the conventional dipole- 
dipole interaction form 



y-dipole 



|R,,fl-3R,,R,, , ^25) 



where R^j = Ri — R^ . 

The commutation relations of the quantum electrody- 
namic fields in the Coulomb gauge are given by 

[E(r), A(r')] = ATTihcA{r - r') (26) 

where the transverse delta function is defined as 



A(r) = y (l - kk) e 
Equivalently 
A(r) = 
so that 



V(27r)3 



kl^l -kk\ ^^k-r / 



(27r)= 



Exphcitly m 
Eqs.(25) and (30) imply 



(27) 



(28) 



(29) 



(30) 



dipole 



~47r/i, • A(R„) • fij for Ry ^ 0. (31) 



The transverse distribution function A(r) removes the 
longitudinal part of a vector field. For example, if 



Af 



P(r) = ^M,'5(r-R,) 



(32) 



represents polarization, i.e. the molecular dipole moment 
per unit volume, then 



P(r) ^ J A(r-r') •P(r')d^r' 

represents the transverse part of the polarization. 
Now, consider the unitary operator 



5 = exp ( ^ / A(r) • P{r)d^r 



(33) 



(34) 



Eq.(26), (33) and (34) imply 

5^E(r)5 = E(r) - 47rP(r) (35) 
so that the radiation Hamiltonian of Eq.(23) reads 

S^HradS = y (|E - 47rP|2 + \B\'')d\. (36) 

Eqs.(17), (23), (31)-(33) and (36) imply that 

S^HradS = Hrad - / E • Pd^ 



dipole 



(37) 



i<i<j<j^ 



It is useful to define the "screened" (or short ranged) 
intermolecular Coulomb potential with the long ranged 
dipole-dipole potential subtracted 



3 



(38) 



Finally, the unitary transformation removes the vec- 
tor potential from the internal degrees of freedom of a 
molecule as in the above Sec. II; i.e. in the dipole inter- 
action limit 



(39) 



The "grand finale" of the algebraic ceremony of this 
section is the assertion that total transformed Hamilto- 



H = S'^HS 
which for J\f molecules takes the form 

7^ = i-/(E2 + B2)rf3r+ 



(40) 



l<i<j<J\f 



(41) 



where the internal molecular Hamiltonian 

'Hmol-j = Hmol-ji^ = 0) + Wmol-j (42) 

is independent of the transverse A. The central Eq.(41) 
of this section follows from Eqs.(17), (24), (37)-(39) and 
(40). 

IV. DICKE-PREPARATA MODEL 

The Hamiltonian in Eq.(41) has four parts: (i) a sum 

of intramolecular Hamiltonians describing the internal ki- 
netic and Coulomb energies of each molecule, (ii) a sum 
of screened (short ranged) intermolecular Coulomb po- 
tentials between neighboring molecules, (iii) the radia- 
tion field energy and (iv) a linear coupling between the 
dipole moment per unit volume and the electric field. 
Within the dipole approximation, Eq.(41) is exact but 
not yet amenable to rigorous mathematical solution. For 
the purpose of obtaining reasonable answers for (say) the 
free energy, further simplifications must be made. 

For the Dicke-Preparata model, the approximations 
are as follows: (i) The s\im of intramolecular Hamilto- 
nians is modeled by a sum of "two-level molecules" , each 
described by Pauli matrices (fx;j, crj^;j, (Tzy); 



(43) 



(ii) The sum of screened intermolecular Coulomb poten- 
tials between neighboring molecules is neglected 



0. 



(44) 



(iii) The radiation field energy is modeled by a single 
photon oscillator mode 



(45) 



(iv) The linear coupling between the dipole moment per 
unit volume and the electric field is modeled by the res- 
onant photon creation and annihilation terms 



/ 



E • Pdrr ^ fl 



,3 ~ / lirtiujo 



n 



^(6a+;, + &ta_.,.), (46) 



where '2a±-j = ax-j ± i^^yj and p. is the excitation ma- 
trix element of the electric dipole operator of a single 
molecule. Putting the model replacement Eqs.(43)-(46) 
into Eq.(41) yields the Dicke-Preparata model 



(47) 



where 



nDP=eS, + n^ob^b + flJ'^^^ {bS++b^S-), (48) 



1 ^ 



and 



1 ^ 



(49) 



(50) 



The Dicke-Preparata model may be written in terms 
of the oscillator coordinate Q and momentum P using 



^ ^ (P - lUgQ) ^ (P + iuj^Q) 



yielding 



where the vector h is given by 



h = M\/ {Pi + woQj) + £k. 



(51) 



(52) 



(53) 



The free energy T of the model is determined by 

Here, the trac-c; over the oscillator degree of freedom may 
be taken as being classical with virtually zero error in the 
limit A/' — > oo; i.e. 



J — oo J — oo 



\ 2nn J 
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Since 



tr e 



-h-S /ksT 



2 cosh 



2kBT 



and since Eq.(53) implies 



it follows from Eqs.(52), (53), and (55)-(57) that 



OO fOO 



— OO 'J — OO 



-Qlk.T ( dPdQ \ 
V 2nh ) 



where 



w — ksT In 2 cosh ■ 



2kBT 

the oscillator energy per molecule is 



w 



2Af 



and the number of molecules per unit volume n is 
n= (Af/n). 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 



In the thermodynamic limit of Eqs.(58)-(61), the min- 
imum free energy principle is that 



lim (T/M) = inf f{w,n,T), 

Af^oc (0<Mi<oo) 

where 

f{w, n,T) = w — ksT In (2 coshx(w, n, T)) 

and 



x{w,n,T) 



2fc^r 



(62) 



(63) 



(64) 



Let 'w{n,T) > represent the mean oscillator energy 
per molecule as a function of the molecular density n 
and the temperature T. If u; > 0, then the system is 
in a super-radiant phase. If = 0, then the system 
is in a normal (incoherent radiation) phase. The free 
energy per molecule f(w,n,T) is at a minimum for w = 
w. Setting {df/dw)nT = yields the following implicit 
function equation for w: 



kBTx{w,n,T) 



tanhx(w, n, T). 



(65) 



The phase diagram in the (n, T) plane follows from 
setting w{n,Tc) = and solving for the critical temper- 
ature Tc{n) as a function of density. The super-radiant 
phase exists for T < Tc{n) and the normal phase exists 



for T > Tc{n). The critical temperature Tc{n < nc) = 
where the critical density 



27r/i2 



Analytically, 

ksTcin > 71c) = 



ln[{n + nc)/{n - Uc)) 
A plot of the phase diagram is shown below. 

Phase Diagram 



(66) 



(67) 
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FIG. 2. Shown is the phase diagram in the (n, T)-plane, 
where is defined in Eq.(66). The two phases are separated 
by the critical temperature curve {kBTcin) / e). 



V. CONCLUSIONS 

The doubts about the physical reality of the super- 
radiant phases were discussed in great detail by I. 
Bialynicki-Birula and K. Rzazewski (23|, who came to 
the remarkable (and dubious) conclusion that the electric 
dipole-photon interaction could be gauged away. This in- 
correct conclusion was based on using a vector potential 
which is independent of both space and time. It is no 
wonder that such a pure "gauge field" had zero effect on 
the thermodynamic equations of state. 

In reality, the transverse electric field in the Coulomb 
gauge should be calculated using either 

1 9A 

E = — , (Heisenberg or Classical Picture), (68) 

c at 



E 



[iJ, A], (Schrodinger Picture). (69) 
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The classical canonical transformation Eq.(12) of I. 
Bialynicki-Birula and K. Rzazewski |^ (in the dipole 
approximation) should have been properly performed as 
in om- Eqs.(5) and (6). The last term on the right hand 
side of both our Eqs.(5) and (6) were incorrectly omit- 
ted in Eq.(12) of I. Bialynicki-Birula and K. Rzazewski 
psf . This incorrect classical canonical transformation 
eliminated the electric field E and thereby the electric 
dipole-photon interaction. 

A proper procedure is to employ a fully quantum elec- 
trodynamic framework to obtain a non-trivial electric 
field E before going to the classical electrodynamic os- 
cillator trace evaluation. From this rigorous Schrodinger 
picture point of view, Eq.(12) of I. Bialynicki-Birula and 
K. Rzazewski [|3| is still incorrect because Eqs.(B2-B3) 
in our Appendix B Eqs.(B2-B3) were not taken into ac- 
count in that work P3| . 

In any or all pictures, the electric dipole-photon inter- 
action cannot be described by gauging away the electric 
field E. When the Hamiltonian is expressed in terms of E 
and B, it is then properly gauge invariant and non-trivial 
changes very surely appear in the thermal equations of 
state. For example, the one-photon loop corrections to 
the equations of state have long been known to produce 
Casimir-van der Waals forces in condensed matter. Such 
forces cannot possibly be gauged away and I. Bialynicki- 
Birula and K. Rzazewski p3[ are thereby in error. A 
further error is found in K. Rzazewski, K. Wodkiewicz 
and W. Zakowicz |^^, where it is asserted that the super- 
radiant phase transition is null and void due to quadratic 
coupling in the photon creation and annihilation opera- 
tors. Had the canonical transformation been carried out 
correctly, as in the present work, then it would have been 
noted that the coupling between the electric dipole mo- 
ment and the photons in the gauge invariant Eq.(41) is 
strictly linear. The super-radiant phase transition re- 
mains intact. 

Preparata has recently proposed several examples of 
super-radiant phases in a recent review |30[| . Water is per- 



haps the most interesting example 31|. In future work 
it would appear interesting to relate the super-radiant 
transition to an instability of the incoherent radiation 
phase which produces the conventional theoretical and 
experimental QED Casimir effects. QED effects at the 
one photon loop level have long proved their importance 
in the understanding of the long ranged intermolecular 
forces |9|. 

MATHEMATICAL APPENDIX A 

Starting from the Schrodinger equation 

one seeks a solution of the form 

\m = u{t) m)) , {A2) 



where W{t) = U-'^{t). Putting Eq.(A2) into Eq.(Al) 
yields 



\m) 



or equivalently 



■H{t) \m) 



(A3) 



(A4) 



where 



■H{t) = uHt)H{t)u{t) - ihU^it)^:^. {A5) 

Eq.(A5) is the required time dependent canonical trans- 
formation. 

MATHEMATICAL APPENDIX B 

From Eqs.(l) and (14) it follows that 

S^PaS - Pa + -A. (Bl) 

c 

Eqs.(Bl) and (11) then imply 

S^Hmoi{A)S^Hmoi{A = 0)- {B2) 

Thus, the vector potential in the dipole interaction theory 
may be "gauged away" . However, the dipole interaction 
returns to the full Hamiltonian since Eqs.(9) and (14) 
imply that 



27r 



and 



2tt 



{B2) 



m 



fihujk 

From Eqs.(7), (B2) and (B3) it follows that 

S^HradS = - M • E + Wrnol (B4) 

where the operator electric field E is given in Eq.(lO), and 
the molecular self energy Wmoi is given in Eq.(17). Thus 
the same canonical transformation which gauges away 
the vector potential in the molecule Hamiltonian, brings 
back the electric dipole interaction — /i • E when applied 
to the radiation Hamiltonian via Eq.(B4). Finally, from 
Eqs.(B2) and (B4) it follows that 



Hrad + HrnoliA = 0) + Wrnol - ' E, 



(55) 
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and Eqs.(16) and (18) then follows from Eq.(B5). The 
proof of Eq.(18) has been completed. 
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